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ABSTRACT 


Weiss  and  Weinstein  have  recently  derived  a  new  lower  bound  on  random 
parameter  estimation  error.  This  report  makes  use  of  this  bound  to  develop  a 
criterion  for  comparing  two  sparse  array  configurations.  The  arrays  are  first  adjusted 
in  size  so  that  their  direction-finding  (DF)  accuracies  at  high  element  signal-to-noise 
ratio  (ESNR)  are  the  same  when  averaged  over  all  possible  directions  of  arrival  in  the 
hemisphere.  Then  their  threshold  ESNRs  are  calculated  using  the  Weiss- Weinstein 
bound.  The  array  having  the  lower  threshold  ESNR  is  judged  the  better  array.  As 
an  example  of  the  application  of  the  bound,  three  7-element  arrays,  a  circle,  a  spiral, 
and  an  L-shaped  array,  are  compared. 

The  bound  is  also  extended  to  the  case  of  switched  subarrays,  where  data  are 
collected  from  several  subarrays  sequentially  in  time  and  processed  as  a  unit.  The 
purpose  is  to  achieve  DF  accuracy  approaching  that  of  the  full  array  while  reducing 
the  number  of  signal  processing  channels.  An  example,  using  five  receivers  on  the 
7-element  spiral  array,  shows  that  this  is  indeed  possible. 
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1.  INTRODUCTION 


The  purpose  of  this  report  is  to  develop  a  method  for  comparing  two  sparse  planar  array 
configurations.  The  criterion  for  comparison  is  based  on  narrowband  single-signal  direction  finding 
(DF)  accuracy  in  the  presence  of  spatially  white  noise  (i.e.,  uncorrelated  from  receiver  to  receiver). 
The  directions  of  arrival  are  assumed  to  be  confined  to  a  he.aisphere,  and  DF  error  is  averaged  over 
all  possible  directions  of  arrival  in  the  hemisphere. 

According  to  sampling  theory,  the  aperture  illumination  must  be  sampled  about  every  half¬ 
wavelength  in  both  spatial  dimensions  in  order  to  preserve  all  the  available  information  about 
the  directions  of  the  signals.  Each  sample  requires  an  antenna  element  and  a  receiver  channel. 
Unfortunately,  system  cost  increases  very  rapidly  with  the  number  of  receiver  channels  used  in  the 
array. 

The  performance  of  an  array  is  characterized  by  its  array  pattern,  which  is  the  response 
of  the  array,  phased  to  point  in  a  particular  direction,  to  a  signal  arriving  from  another  direction. 
When  the  array  satisfies  the  half-wavelength  sampling  criterion,  the  array  pattern  has  a  well-defined 
main  beam  and  uniformly  low  sidelobes.  As  elements  are  removed  from  the  array,  the  sidelobes 
grow  in  size  and  number.  A  high  sidelobe  means  that  there  are  two  directions  having  very  similar 
array  response  vectors.  In  the  presence  of  noise  and  calibration  uncertainties,  confusion  between 
the  mainlobe  and  the  sidelobes  can  result.  If  the  number  of  signals  is  known,  or  can  be  reliably 
estimated,  this  confusion  manifests  itself  as  occasional  large  DF  errors;  if  not,  it  causes  false  alarms 
as  well. 

The  value  of  element  signal-to-noise  ratio  (ESNR)  at  which  these  large  errors  begin  to  occur 
with  significant  probability  is  called  the  threshold  ESNR  of  the  array.  If  two  arrays  have  equal 
accuracy  at  high  ESNR,  the  array  having  the  lower  threshold  ESNR  is  deemed  the  better  array. 

The  Weiss- Weinstein  (WW)  bound  is  a  lower  bound  on  the  mean  squared  estimation  error 
for  random  parameters  with  known  a  priori  distributions.  As  will  be  shown,  it  predicts  remarkably 
well  the  threshold  ESNR  of  the  maximum  likelihood  direction  estimator.  This  quantity  would 
otherwise  have  to  be  determined  by  simulation. 

Equations  for  the  bound  are  developed  for  the  planar  array  problem,  and  various  issues  that 
arise  in  its  calculation  are  addressed.  A  scalar  mecisureof  accuracy  is  proposed  for  two-dimensional 
DF  problems,  and  the  Cramer-Rao  (CR)  factor  oi  an  array  is  defined.  Arrays  having  the  same  CR 
factor  have  the  same  DF  accuracy  bound  at  high  ESNR.  As  an  example  of  the  application  of  the 
bound,  three  7-element  arrays,  a  circle,  a  spiral,  and  an  L-shaped  array,  are  compared. 

The  bound  is  extended  to  the  case  of  switched  subarrays,  where  data  is  collected  from  several 
subarrays  sequentially  in  time  and  processed  as  a  unit.  The  purpose  is  to  achieve  DF  accuracy 
approaching  that  of  the  full  array  while  reducing  the  number  of  signal  processing  channels.  An 
example,  using  five  receivers  on  the  7-element  circular  array,  shows  that  this  is  indeed  possible. 
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The  main  advantage  of  adaptive  arrays  is  their  ability  to  adapt  their  array  patterns  to  the 
signal  environment,  amplifying  the  signals  of  interest  while  simultaneously  redui  ing  other  unwanted 
signals.  Nevertheless,  as  a  first  step  toward  a  rationale  for  adaptive  array  comparison,  arrays  will 
be  compared  on  the  basis  of  their  single-signal  DF  accuracy,  assuming  ideal  (i.e..  omnidirectional) 
antenna  element  patterns.  If  this  approrch  provides  some  insight  into  the  problem,  then  perhaps 
it  can  be  extended  to  more  complicated  scenarios. 

The  following  notational  conventions  are  used.  Bold  lowercase  letters  represent  (column) 
vectors,  while  bold  uppercase  letters  represent  matrices.  No  particular  distinctloi\  is  made  between 
real  and  complex  quantities;  the  context  should  make  it  clear.  The  transpose  of  a  matrix  A  is  A^: 
its  complex  conjugate  is  A*;  its  conjugate,  or  Hermitian,  transpose  is  A^.  The  determinant  of  A 
is  denoted  by  det  A.  I  denotes  the  identity  matrix;  1  denotes  a  vector  having  all  elements  ecjual 
to  1.  The  notation  [x]  denotes  the  diagonal  matrix,  the  nonzero  elements  of  which  are  those  of  the 
vector  X.  The  inequality  A  >  B  means  that  the  matrix  A  -  B  is  nonnegative  definite.  .All  matrix 
inverses  (including  scalars)  that  appear  are  a,ssumed  to  exist,  unless  otherwise  noted. 


2.  THE  PARAMETER  ESTIMATION  PROBLEM 

DF  is  an  example  of  a  parameter  estimation  problem.  An  excellent  discussion  of  such  problems 
is  available  in  Van  Trees  [1].  The  following  brief  description  is  given  primarily  to  introduce  some 
terminology  and  notation. 

Observations  of  a  (scalar-  or  vector- valued )  random  process  are  made  in  order  to  determine 
which  of  a  set  of  possible  probability  density  functions  (pdfs)  generated  them,  if  the  set  of  possible 
pdfs  is  finite  (usually  sniall),  the  problem  is  called  a  hypothesis  testing  problem;  otherwise,  it 
is  called  a  parameter  estimation  problem.  Parameters  are  quantities  that  may  be  thought  of  as 
indexing  the  set  of  possible  pdfs,  i.e.,  knowing  the  values  of  the  parameters  is  equivalent  to  knowing 
which  pdf  is  true.  The  vector,  the  elements  of  which  are  the  parameter  values,  is  railed  the  pararmttr 
vector. 

To  best  model  the  physical  situation  at  hand,  the  elements  of  the  parameter  vector  may  be 
treated  either  as  unknown  constants  or  as  random  variables  with  a  joint  pdf  of  their  own.  When 
they  are  treated  as  unknown  constants,  the  set  of  possible  pdfs  is  usually  written  as  pz(z;  u),  where 
z  denotes  the  observation  and  u  denotes  the  parameter.  When  the  parameters  are  treated  as 
random,  the  jotnt  pdf  for  z  and  u  is  written  as  pz,u(z.u),  the  a  priori  pdf  for  u  is  the  marginal 
density 


Pu(u)  =  j  Pz,u(2.u)dz, 

and  similarly  for  z,  and  the  a  posteriori  density  for  u,  i.e..  the  pdf  for  u  given  the  values  of  the 
observations,  is  the  conditional  density 

,  ,  ,  Pz.ul^-tt)  Pz|u(2|u)pu(u) 

Pu|z(u|z)  =  =  - -7-T - • 

Pz(z)  Pzl*) 

A  similar  notation  is  used  for  expectations  of  random  quantities,  e.g.,  the  conditional  expectation 
of  a  function  /(u)  of  u,  given  the  value  taken  by  z,  is  written 

Jpu|z{f(u)}  =  j  f(u)p(ulz)du. 
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Sometimes,  when  it  is  clear  from  the  context  what  is  intended,  the  subscript  on  p  will  be 
omitted,  as  in  the  previous  equation.  In  this  report,  the  parameter  vector  will  be  treated  as 
random.  An  estimate  of  u  is  a  function  only  of  the  observations,  u(z).  The  estimation  error. 
u{z)  -  u,  is  a  random  vector,  and  so  is  completely  characterized  by  the  join;  pdf  pi  z  u).  In  many 
cases,  it  ran  be  argued  *hat  the  error  vector  is  approximately  Gaussian  and  so  is  well  characterized 
by  its  mean  value  and  covariance  matrix,  or  by  its  mean  squared  error  matrix 

i’2.u{(u(z)  -  u)(u(z)  -  u)^}. 
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3.  THE  WW  BOUND 


Weiss  and  Weinstein  [2,3,4]  have  recently  derived  a  bound  on  the  mean  squared  error  '  in 
estimating  random  parameters.  The  bound  states  that  [4,  Equation  (18)] 


£z,u(u-u)(u-u)'^  >Hr-'H^  =  n,  (1) 

where  is  the  vector  estimate,  is  the  true  parameter  vector,  is  a  matrix,  the 

columns  h,  of  which  represent  perturbations  in  the  true  parameter  vector  u,  and  is  a  matrix 

to  be  defined  below  in  terms  of  the  likelihood  ratio  of  the  observations  when  the  true  parameter 
vector  is  u  +  h  and  u,  respectively. 


L(z;  u  +  h,  u)  = 


p(2,u  +  h) 
p(z,u) 


(2) 


The  likelihood  function  is  defined  only  on  the  set  U  of  values  of  u  for  which  p(z,  u)  >  0.  Define  a 
real  random  vector  w,  the  ith  element  of  which  is  given  by 


w, 


A  L*/^(z;  u  +  hj,u)  -  u  -  hi,u) 

£{L‘/2(z;u  +  hj.u)} 

The  matrix  F  is  then  simply  the  covariance  matrix  of  w, 

r  =  £{ww^}. 


(3) 


It  is  clear  from  this  representation  that  F  is  nonnegative  definite;  hereafter  it  is  assumed  to  be 
positive  definite. 

The  number  T  and  nature  of  the  perturbation  vectors  in  H  are  completely  arbitrary.  Every 
choice  of  them  produces  a  lower  bound  on  the  mean-squared  error  matrix.  The  addition  of  another 
perturbation  vector  can  only  increase  the  bound,  in  the  sense  that  if  the  bounding  matrix  is  fl 
before  and  il'  after  a  perturbation  vector  is  added,  then  SI'  —  SI  >  0.  A  proof  of  this  result  is  given 
in  Appendix  A. 


'Their  bound  is  actually  more  general  than  the  one  described  here.  See  Weiss  and  Weinstein  [4] 
for  the  more  general  form. 
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4.  THE  SINGLE-SOURCE,  IDEAL-ARRAY  DF  PROBLEM 


The  mathematical  model  of  the  array  DF  problem  for  narrowband  sources  that  will  be  used 
here  has  been  described  in  many  papers  (see,  e.g.,  Johnson  [5]).  A  brief  description  will  be  given 
here  in  order  to  introduce  the  notation  to  be  used.  Attention  is  restricted  here  to  the  case  of  a 
single  narrowband  source,  so  that  the  complex  envelopes  of  the  signals  received  at  each  of  the  M 
array  elements  at  time  t  form  a  complex  M  vector 

z(t)  = /3(t)v(u)  +  n(t),  (4) 

where  u  is  a  real  unit  vector  pointing  from  the  array  origin  to  the  source,  i3(t)  is  the  complex  i 

amplitude  of  the  incident  wavefront  at  the  array  origin.  v(u)  is  a  complex  M  vector,  the  array 

response  vector  for  direction  u,  containing  the  complex  signal  components  received  at  each  of  the 

elements  divided  by  j3(t),  and  n(f )  is  a  complex  M  vector  containing  the  noise  components  at  each 

element.  The  data  consists  of  samples  of  z(t)  taken  at  times  (ti,t2i  •  •  •  <  <^’)•  For  brevity,  all  vector 

values  at  time  t„  will  be  indexed  simply  by  n.  It  is  assumed  that  the  time  between  samples  is  large  I 

enough  so  that  the  noise  vectors  n„  and  also  the  complex  signal  amplitudes  are  statistically 

independent.  The  noise  vectors  are  modeled  as  circular  complex  Gaussian  random  vectors  with 

zero  mean  and  covariance  matrix  1.^  The  complex  signal  ampbtude  /3„  is  also  modeled  as  a  circular 

complex  Gaussian  random  variable  with  zero  mean  and  known,  constant  variance  F.^  With  these 

assumptions,  the  data  vectors  z„  are  independent,  circular  complex  Gaussian  random  vectors  with  ^ 

zero  mean  and  common  covariance  matrix  C  =  Pv(u)v^(u)  +  I. 

The  signal  source  is  assumed  to  be  far  away  from  the  array,  so  that  the  incident  wavefront  is 
essentially  a  plane  wave.  If  the  mth  ideal  element"*  is  located  at  position  dn,  relative  to  the  array 
origin,  its  array  response  is  expfjfcu^d^),  where  k  =  27r/A  is  the  wave  number.  The  array 

response  vector  (transposed)  is  I 

v^(u)  =  (exp(jA:u^di),exp(j7:u^d2),...,exp(;fcu^dM))-  (5) 


^All  2M  real  and  imaginary  parts  of  the  entries  of  nr.  are  statistically  independent  and  have  equal 
variances  of  1/2. 

^If  the  variance  P  is  unknown,  it  becomes  an  additional  parameter  to  be  estimated,  and  must  be 
assigned  an  a  priori  probability  density  if  one  wishes  to  use  the  WW  bound. 

"‘An  ideal  element  measures  the  complex  amplitude  of  the  incident  wavefront  with  an  ideal  omni-  • 

directional  pattern. 
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Note  that  ||v(u)||^  =  M.  Because  the  signal  magnitude  at  each  element  is  \(i\  and  the  noise  variance 
at  each  element  is  unity,  P  represents  the  ESNR  at  each  element. 

The  array  is  assumed  to  be  planar.  This  assumption  affects  the  analysis  in  that  it  simplifies 
the  form  of  the  derivative  with  respect  to  u  of  v(u).  This  means  that  the  z  component  of  all  the 
vectors  are  zero,  and  so  they  and  the  vector  u  may  be  regarded  as  two-dimensional.  It  is  also 
assumed  that  the  direction  of  arrival  is  confined  to  the  hemisphere  above  the  plane  of  the  array. 
This  eliminates  any  ambiguity  in  the  z  component  of  the  vector  u. 

There  are  thus  two  real  parameters  to  be  estimated,  the  x  and  y  components  of  u.  An  a 
priori  pdf  must  therefore  be  chosen  for  these  variables.  To  keep  the  computation  of  the  bound  as 
simple  as  possible,  this  function  is  chosen  to  be  constant  on  the  unit  disk,  i.e., 


P{U)  =  p(Ui,Uy)  = 


1/rr  if  -1-  <  1 

0  otherwise. 


This  density  implies  that  the  marginal  probability  density  of  signals  in  elevation  angle  is  sin  2<^. 
If  the  signal  directions  were  uniformly  distributed  in  solid  angle  over  a  hemisphere,  the  resulting 
density  in  elevation  angle  would  be  cos<p,  producing  more  signals  near  the  horizon  and  less  at  higher 
elevation  angles.  In  many  applications,  this  would  be  a  more  realistic  model;  however,  it  would 
greatly  complicate  the  computation  of  the  bound. 


•  • 
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5.  CALCULATION  OF  THE  WW  BOUND  ^ 

The  first  step  in  computing  the  bound  is  to  write  down  the  joint  probability  density  of  the 
measured  data  and  the  parameters,  which  is  the  product  of  the  a  priori  density  on  u  and  the 
conditional  density  of  the  data  when  u  is  known.  This  conditional  density  is  circular  complex 
Gaussian,  with  covariance  matrix  § 


C(u)  =  F{zz"}  =  I  +  Pv{u)v"(u). 


(7) 


The  inverse  of  C  is 


C-‘(u)  =  I- 


v(u)v"(u), 


MP+l 


(8) 


and  its  determinant  is 

detC(u)  =  MP+l. 


The  joint  probability  density  of  the  data  z  =  {zn,n  =  1,  and  the  parameters  is 

N 

p(z,u)  =  p(z|u)p(u)  =  lir^(MP  +•  Ijj'^^expi-  52znC"’(u)zn}p(u).  (9) 

n=l 


The  likelihood  ratios  required  in  the  evaluation  of  the  bound  are 


r,  ..  ,  ^  P(Z.U+h,) 

Lz;u+h.,u)  =  ^ 

p  z,u 


P(u  +  hi) 

p(u) 


exp 


MP+l 


^(|z"v(u  +  h,)|2  -  lz^v(u)|2) 


Note  that  these  likelihood  ratios  are  defined  only  on  the  set  U  (the  unit  disk)  where  p(u)  >  0. 

All  that  is  needed  to  compute  the  F  matrix  [Equation  (3)]  required  for  the  WW  bound  is  the 
evaluation  of  two  expectations.  The  first  is 
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h,,u)}  =  ^  dup(u) 

•  J  ‘^exp|^-jj^^^^{|z"v(u+h,)|^  -  |z"v(u)|^)| 
[ff'^(MP  + ir'^expj- x;;z«c-*(u)z„| . 

The  Itist  two  factors  under  the  z  integral  above  are  of  course  just  p(z|u).  The  right  side  simplifies 


^dup'/^(u  + h,)p‘^^(u)  j  dz[n^{MP  +  1)]  ^  exp{- z"Qz„}, 


where 


Q  =  I  -  ^-^^^^(v(u  +  h.)v"(u  +  h,)  + v(u)v"(u)).  (11) 

The  z  integration  reduces  to  a  product  of  N  integrals  of  the  form 

[i:^{MP  +  1)]"*  j  (iz„exp{-z^Qz„}  =  ((A/P  +  l)detQ]"'. 

It  is  shown  in  Appendix  C  that 

(A/P+l)detQ  =  i  +  _J^|!_(l_|v«(u)v(u  +  h.)|VM2) 

=  1  +  -  |v"(0)v(h,)lVM^).  (12) 

where  v(0)  =  1.  Note  that  detQ  is  not  a  function  of  u,  and  so  the  two  integrals  in  Equation  (10) 
are  separable.  The  final  result  is 

£{T’'^^(z;u  +  hj,u)}  =  [(A/P  +  IjdetQ]''^  /  du  p^^^(u  +  h,)p’^^(u).  (13) 

J  u 


The  second  required  expectation  is 

u  +  h,,  u)£*/^(z;  u  +  hj,  u)}  = 

■lTr^(MP+\)]~^  f  dz  exp|-^^z"C-’lu)z„| 

expj^^-^—  5^(|z"v(u  +  hJpH-  |z"v(u  +  hj)p  -  2|z"v(u)|^)| , 

which  reduces  to  [compare  Equation  (10)] 

j^dup’^^(u  +  hi)p*''^(u  +  hj) y^z  [x*^(A/P+ l)]~^exp{- 5^^z"Rz„},  (14) 

where  [compare  Equation  (11)] 


R  =  I  -  ^-;jj^;^(v(u  +  h,)v"(u  +  h,)  + v(u  +  h^)v"(u+  hj)). 

The  same  sequence  of  steps  that  produced  Equation  (13)  now  leads  to 
E{  L’/^(z;  u  +  h, ,  u)L‘/*(z;  u  +  hj,  u)}  = 

l(MP  +  l)det  R]~^  du  p‘^^(u  +  hj)p‘^^(u  +  h^) 
where  [compare  Equation  (12)] 

(MP+UdetR  =  l  +  ^^^^^|;^(l-lv"(u  +  hi)v(u  +  h;)|VA/*) 


=  1  + 


J^(.-|v''(h.Wh,)|VM=). 


(15) 


(16) 


(17) 


Again  detR  is  not  a  function  of  u,  and  so  the  two  integrals  in  Equation  (14)  are  separable,  yielding 
Equation  (16).  Note  that  the  determinant  detQ  given  by  Ek}uation  (12)  is  just  a  special  case  of 
Equation  (17). 
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In  order  to  obtain  a  compact  expression  for  the  F  matrix  that  defines  the  WW  bound,  it  is 
convenient  to  define  the  functions 


AN(h^,h,)  =  [(MP+l)detR]-'^, 

(18) 

B(h,,hj)  =  ^dup’/^(u+ h,)p‘/^(u  + hj). 

(19) 

and  their  product 

C(h,,hj)  =  i4A/(h,,hj)5(h,,hj). 

Note  that  only  An  depends  on  the  signal-to-noise  ratio  P  and  the  number  of  samples  N.  Note 
also  that  because  the  domain  of  integration  in  Equation  (19)  is  the  set  U,  the  region  of  support  in 
Equation  (19)  is  the  intersection  of  three  unit  disks,  not  two  (see  Appendix  D). 

In  terms  of  these  fui.ctions,  the  ijth  element  of  the  F  matrix  of  Equation  (3)  takes  the  form 

r  -  C(h.'  h.)  +  C(-h„  -h,)  -  C(h.,-h,)  -  C(-h.,h,) 

C(h.,0)C(h„0) 


It  is  clear  from  their  definitions  that  An  and  B  have  the  common  property 
AAf(hi,hj)  =  A/v(hj,hi), 

B(hi,hj)  =  Bih„h,), 

and  that  A/v  is  an  even  function  in  the  sense  that 
AA/(-hj,  -hj)  =  A/v(h,,  hj). 

If,  as  has  been  assumed  here,  the  parameter  set  U  has  the  property  that  if  u  is  in  U,  — u  is  also  in 
{/,  and  that  the  a  priori  density  p(u)  is  an  even  function,  i.e.,  p(  — u)  =  p(u),  then  Equation  (20) 
simplifies  to 


2[C(h.,h,)-C(h.,-h,)] 

C(h.,0)C(hj,0) 
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6.  A  SCALAR  MEASURE  OF  DF  ACCURACY 


A  planar  array  can  estimate  two  real  angular  parameters  (here  Uj.,Uy)  that  define  a  unique 
direction  in  a  three-dimensional  hemisphere.  The  DF  error  is  thus  a  two-vector  having  a  two- 
dimensional  probability  density.  Different  array  geometries  will  produce  distributions  having  dif¬ 
ferent  shapes.  In  order  to  compare  two  such  arrays,  their  distributions  must  be  characterized  by  a 
single  real  number. 

A  convenient  choice  is  the  area  of  the  ellipse  of  concentration  of  the  error  distribution  (e.g., 
see  Cramer  [6,  p.283]).  A  probability  density  for  the  error  uniformly  distributed  over  the  ellipse  of 
concentration  would  have  the  same  mean  squared  error  matrix  as  the  actual  a  posteriori  density. 

This  area  of  the  ellipse  of  concentration  is  given  by  4irv'QetS,  where  S  is  the  mean  squared  error  ^ 

matrix  of  the  distribution.  The  array  having  the  smallest  area  of  concentration  is  deemed  the  best. 

It  is  convenient  to  plot  this  error  area  on  a  logarithmic  scale. 

In  the  limit  as  ESNR  goes  to  zero,  the  observations  provide  no  information  and  the  mean 
squared  error  matrix  becomes  that  of  the  a  priori  density.  This  has  S  =  1/4,  so  that  det  ii  =  1/16, 
and  its  ellipse  of  concentration  has  area  ir.  This  is  the  value  the  scalar  measure  should  approach  I 

for  low  ESNR. 


» 


» 


( 
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7.  THE  CR  FACTOR  OF  AN  ARRAY 

General  expressions  for  the  Fisher  information  matrix  F  for  nonrandom  parameter  estimation 
are  given  in  DeLong  [7]  for  the  case  of  a  multiple  signal  scenario  in  which  each  signal  has  a 
single  unknown  directional  parameter.  Both  deterministic  and  random  Gaussian  signal  models  are 
considered.  In  the  deterministic  model,  the  complex  amplitudes  may  be  either  known  or  unknown. 
In  the  Gaussian  model,  the  covariance  matrix  of  the  signal  complex  amplitudes  may  be  either 
known,  or  unknown  but  diagonal  (uncorrelated  signals).  The  Gaussian,  known  covariance  bounds 
derived  from  this  information  matrix  (assuming  unbiased  estimates)  exhibit  unusual  behavior  in  the 
multiple  signal  case,  which  has  been  verified  by  simulation  of  estimation  algorithms  in  White  [8]. 

The  Fisher  matrix  of  DeLong  [7,  Equation  (2.242)],  for  the  Gaussian-signal,  known-covariance 
problem,  is  extended  to  signals  with  an  arbitrary  number  of  directional  parameters  (e.g.,  the  planar 
array  problem)  in  Appendix  F.  When  specialized  to  the  case  of  a  planar  array  of  omnidirectional 
elements,  it  becomes 


F  =  2k^N 


l  +  MP 


D^(I-±11^)D 


where 


D^  = 


X\  X2 

yi  yz 


is  a  matrix  containing  the  x,y  positions  of  the  elements  of  the  array  in  some  arbitrary  coordinate 
system.  Note  that  F  is  not  a  function  of  u.  This  fact  will  be  quite  useful  below.  For  Gaussian 
signals,  the  quantity  F“’  is  a  lower  bound  on  the  covariance  matrix  of  any  unbiased  estimate  of  u. 
It  is  referred  to  hereafter  as  the  unknown  parameter  CR  bound. 

For  large  signal-to-noise  ratios,  the  regime  of  interest  below.  Equation  (21)  becomes 


F  =  2k^NP 


■)d]^ 


2k^NPA. 


The  array  configuration  appears  in  this  expression  entirely  by  way  of  the  matrix  A.  This  matrix 
is  simpler  than  it  appears  at  first  glance.  If  the  origin  of  the  array  coordinate  system  is  moved  to 
the  center  of  mass  of  the  element  positions  (so  that  D^l  =  0  ),  it  becomes  simply 


A  =  D^D  =  ^ 

m=l 


In  the  limit  of  high  ESNR,  it  has  been  found  empirically  that  the  WW  lower  bound  on 
the  mean  squared  error  matrix,  for  the  problem  of  interest  here,  becomes  essentially  equal  to  the  • 

unknown  parameter  CR  bound  when  the  test  points  are  chosen  properly  (see  Section  8).  That  is, 
at  high  ESNR,  both  bounds  say  that 

S>F'‘.  (23) 

The  area  of  the  ellipse  of  concentration  [6]  associated  with  the  error  matrix  S  is  dTrvdetS.  Because 
det  3  >  (det  F)“\  a  lower  bound  on  this  area  is  therefore 

dx-N/det  S  >  (2iryV/>)-'A2(det  (24)  ^ 

Thus,  for  high  ESNR,  both  the  WW  and  CR  approaches  yield  the  DF  error  bound  given  by 

Equation  (24).  The  CR  factor  o{  the  array  is  defined  to  be  the  quantity  (det  A)'/^,  which  has  the 

dimension  of  length.  Arrays  having  equal  CR  factors  and  operating  at  the  same  frequency  have 

equal  DF  accuracy  bounds  at  high  signal-to-noise  ratios.  It  is  sometimes  convenient  to  include  the  ^ 

factor  in  the  definition  of  CR  factor.  When  this  is  done,  the  CR  factor  is  said  to  have  units 

of  wavelength.  The  caveat  about  operating  frequency  can  be  dropped  from  the  above  statement  if 

the  CR  factors  of  the  arrays  are  equal  in  wavelengths. 

A  similar  quantity  for  line  arrays,  L  =  []gEm=i  developed  in  Evans  et  al.  [9], 

where  it  was  called  root  mean  square  (RMS)  array  length.  However,  this  parameter  has  the  unde-  • 

sirable  property  that  removing  an  element  can  cause  it  to  increase.  It  is  easy  to  show  that  the  CR 
factor  defined  here  must  decrease  when  elements  are  removed. 

Although  the  CR  factor  has  the  dimension  of  length,  an  array  with  constant  CR  factor  must 
shrink  as  elements  are  added.  For  this  reason,  the  terms  length  and  aperture  are  not  used  to  describe 
it.  * 


» 
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8.  CHOICE  OF  TEST  POINTS 


Intuition  suggests  and  simulation  confirms  that  the  abrupt  increase  of  the  estimation  error  as 
signal-to- noise  ratio  decreases  occurs  because  noise  causes  a  spurious  sidelobe  peak  to  exceed  the 
true  peak  of  the  likelihood  function  (or  whatever  function  is  being  maximized  by  the  estimator). 
The  unknown  parameter  CR  bound  takes  into  account  variations  of  the  pdf  only  in  the  immediate 
neighborhood  of  the  true  parameter  value  [see  Equation  (F.l)].  Thus  the  performance  threshold  is 
not  well- predicted  by  this  bound.  As  will  be  seen  below,  the  test  points  of  the  WW  bound  can  be 
chosen  so  as  to  predict  the  performance  threshold  quite  well. 

Another  well  known  bound,  the  Barankin  bound,  is  very  similar  in  form  to  the  WW  bound 
and  also  involves  the  use  of  test  points.  McAulay  and  Seidman  [10]  applied  the  Barankin  bound 
to  a  pulse-position  modulation  communications  problem  and  suggested  choosing  the  test  points  to 
be  the  locations  of  the  maximum  sidelobes  of  the  objective  function  (the  matched  filter  output). 

It  is  clear  that  the  height  of  a  sidelobe  is  not  the  only  factor  affecting  the  bound.  If  the 
sidelobe  is  very  distant  in  arrival  angle  from  the  mainlobe,  the  probability  that  the  true  direction 
of  arrival  is  such  that  the  sidelobe  is  within  “visible  space,”  i.e.,  that  it  represents  another  valid 
direction  of  arrival,  is  very  small.  Thus  such  a  sidelobe  contributes  relatively  little  to  the  overall 
DF  error.  Sidelobes  near  the  mainlobe  have  a  greater  probability  of  causing  an  error,  but  the  size 
of  the  error  is  smaller. 

Weiss  and  Weinstein  [2,3,4]  assert  that  if  only  incremental  test  points 


are  used,  the  WW  bound  becomes  the  random  parameter  CR  bound  derived  by  Van  Trees  [1,  pp. 72,84] 
as  f  — ►  0.  It  is  shown  in  Appendix  B,  for  a  simple  one-parameter  problem,  that  this  is  not  the  case 
when  the  joint  pdf  p(z,u)  htis  compact  support  in  u  and  is  not  zero  at  the  end  points.  The  WW 
bound  instead  approaches  0  for  this  case.  This  behavior  is  also  exhibited  by  the  computer  program 
used  to  compute  the  bound  for  the  planar  array  problem  of  interest  here. 

Despite  this  fax;t,  it  is  found  empirically  that  there  exists,  for  every  instance  so  far  examined, 
an  optimal  value  of  f  which  yields  the  tightest  bound.  Furthermore,  the  bound  so  obtained  is  ex¬ 
tremely  close  to  that  obtained  from  the  unknown  parameter  CR  bound.  The  empirically-determined 
equation  for  f  is 

lOlogjoC  =  min[0,0.5862  +  (40-  101og,o(Af Mi^))  -  ESNR(dB)  -  30], 
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4, 


* 


where  N,M  and  L  are  the  number  of  looks,  number  of  elements,  and  CR  factor  in  wavelengths, 
respectively,  of  the  array.  Use  of  this  formula  gives  a  two-point  \VW  bound,  whicii  follows  *he  • 

unknown  parameter  CR  bound  closely  for  ESNRs  ranging  from  —10  to  40  dB 

Experience  shows  that  when  the  test  points  of  the  W\V  bound  are  placed  at  the  locations 
of  the  highest  sidelobes  of  the  array  pattern,  only  a  few  test  points  are  needed.  The  addition  of 
more  test  points  has  very  little  effect  on  the  bound.  An  example  of  this  is  shown  in  Figure  1  for  a 


221432-1 


Figure  I  Effect  of  number  of  test  points  on  the  WW  bound. 


7-element  spiral  array  having  a  growth  factor  of  1.1  and  a  CR  factor  of  8  wavelengths,  processing 
24  looks  of  array  data.  (The  spiral  array  geometry  and  parameters  are  described  more  fully  in 
Section  11.)  The  incremental  test  points  are  the  first  two  used;  they  produce  a  smooth  curve 
without  a  break  point,  essentially  the  unknown  parameter  CR  bound.  Addition  of  two  more  test 
points  at  the  two  highest  sidelobe  locations  produces  the  threshold  behavior,  and  adding  more  test 
points  merely  raises  the  bound  slightly  in  the  uninteresting  ESNR  region  below  threshold.  At  high 
ESNR,  the  WW  bounds  agree  with  the  unknown  parameter  CR  bound,  as  discussed  in  Section  7, 
regardless  of  the  number  of  test  points  (>  2)  used. 

The  next  section  discusses  the  array  pattern  and  the  selection  of  sidelobes  in  more  detail. 
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9.  COMPUTING  THE  ARRAY  POWER  PATTERN 

The  array  power  pattern  function  of  an  array  is  defined  to  be  the  squared  magnitude  of  its 
response,  when  steered  in  direction  Uj,  to  a  signal  incident  from  direction  U2,  i  e., 


1 


(/(Ui,U2)  =  |-j^X(Ul)"x(U2)|^ 


(25) 


where  x(u)  is  the  array  response  vector  for  direction  u.  When  the  array  is  composed  of  ideal, 
omnidirectional  elements,  the  response  vectors  x(u)  become  the  vectors  v(u)  given  by  Equation  (5). 
It  is  then  true  that 


5(ui,U2)  =  ff(0,U2  -  ui)  =  p(h),  (26) 

where  h  =  U2  —  Uj.  An  array  is  characterized  by  computing  g{u\,U2)  over  all  possible  pairs  of 
unit  vectors  in  a  hemisphere.  From  Equation  (26),  it  is  seen  that  this  is  equivalent  to  steering  the 
array  to  zenith  (u  =  0)  and  computing  over  all  possible  differences  of  pairs  of  unit  vectors  in  a 
hemisphere.  It  is  clear  that  the  maximum  length  of  the  difference  of  two  unit  vectors  is  2. 

The  global  maximum  of  the  array  pattern  occurs  at  h  =  0.  Secondary  maxima  of  the  pattern 
(which  may  sometimes  be  equally  large)  are  called  the  sidelobes  of  the  array  pattern.  The  location 
of  a  sidelobe  is  specified  by  a  vector  h.  For  a  general  array,  h  is  three  dimensional,  but  for  planar 
arrays  (the  case  considered  here)  its  component  normal  to  the  plane  is  of  no  consequence.  The 
array  pattern  is  then  a  function  of  two  real  variables,  the  i  and  y  components  of  h. 

A  list  of  sidelobe  values  arranged  in  descending  order,  accompanied  by  their  locations  (h  val¬ 
ues),  defines  the  array  as  far  as  the  WW  bound  is  concerned.  When  the  array  is  steered  to  zenith, 
the  h  values  define  the  sidelobe  locations  directly;  otherwise,  they  define  the  difference  vector 
between  the  steering  direction  and  the  sidelobe  direction. 

While  the  array  pattern  is  translation-invariant,  the  effect  of  a  sidelobe  is  not  translation- 
invariant.  When  the  array  is  steered  to  zenith,  the  only  sidelobes  that  are  important  are  those  lying 
within  the  unit  circle  in  the  h  plane.  This  is  sometimes  called  jisible  space,  because  it  contains  all 
the  physically  possible  directions  of  arrival.  When  the  array  is  steered  to  direction  Uo,  the  main 
peak  moves  there,  and  the  sidelobes  move  with  it.  However,  visible  space  is  still  the  unit  circle. 
Thus,  as  the  beam  is  steered  in  different  directions,  some  sidelobes  vanish  from  visible  space,  while 
others  enter.  DF  accuracy  is  thus  very  much  a  function  of  steering  direction. 

The  WW  bound  overcomes  this  complexity  by  averaging  over  steering  direction,  using  an 
assumed  a  priori  probability  density  for  directions  of  arrival,  as  described  in  Sectioa  4.  A  fixed  set 
of  test  points  is  used;  whether  or  not  they  make  a  contribution  to  the  bound  is  determined  by  the 
factors  5(h,,hj)  defined  by  Equation  (19). 
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When  computing  the  array  pattern,  one  must  choose  a  wavelength.  This  should  be  the 
smallest  wavelength  of  interest.  As  the  wavelength  is  increased,  the  sidelobes  of  the  arr  ,y  pattern 
move  radially  outward  in  direct  proportion.  The  sidelobe  Ust  obtained  for  a  chosen  wavelength  can 
thus  be  used  as  a  source  of  test  points  for  the  WVV’  bound  computation  at  any  longer  wavelength. 


10.  VALIDATION  OF  THE  WW  BOUND 


A  single-signal  DF  simulation  was  carried  out  in  order  to  evaluate  the  tightness  of  the  bound. 
The  DF  algorithm  used  in  the  simulation  was  the  maximum  likelihood  estimator,  discussed  in 
Appendix  E.  The  array  configuration  was  a  5-element  spiral^  with  a  growth  factor  of  2'^^  and  a 
CR  factor  of  4  wavelengths. 

Each  simulation  point  is  the  result  of  averaging  1000  Monte  Carlo  trials.  On  each  trial,  a 
point  (u^.Uy)  was  selected  at  random  from  the  unit  disk,  representing  the  true  direction  of  arrival. 
Using  these  values  in  the  data  model  of  Equation  (4),  N  =  24  data  vectors  (looks)  were  generated 
and  used  to  compute  the  quantity 


n=l 


which  is  monotonically  related  to  the  likelihood  function  (see  Appendix  E).  The  maximum  value 
of  this  function  on  the  unit  disk  |u|  <  1  was  found  by  a  two-dimensional  search. 

The  results  are  shown  in  Figure  2.  Note  that  two  results  are  given  for  ESNRs  of  15  and 
17  dB.  This  is  to  emphasize  the  fact  that,  near  threshold,  the  increase  in  the  error  is  due  to  a  few 
occurrences  of  a  rare  event  (confusion  of  the  main  response  with  a  sidelobe  response).  In  the  first 
se»  of  1000  trials,  this  event  did  not  occur.  A  closer  match  would  presumably  have  been  obtained 
by  conducting  more  trials,  but  the  intensive  computation  required  by  the  two-dimensional  search 
makes  this  impractical.  The  implication  is  that  simulation  is  not  an  efficient  way  to  determine 
threshold  ESNR.  In  view  of  the  difficulties  inherent  in  such  a  simulation,  the  bound  may  be  said 
to  accurately  predict  the  threshold  ESNR  of  this  array. 


^For  a  more  complete  description  of  the  spiral  configuration,  see  Section  11. 


Figure  2  Comparison  of  stmulaiion  results  with  WW  hound. 


11.  COMPAMSON  OF  THREE  ARRAY  GEOMETRIES 


This  section  shows  how  the  WW  bound  can  be  used  to  compare  the  DF  accuracy  of  different 
array  configurations.  The  three  arrays  to  be  studied  have  the  same  CR  factor  and  therefore  the 
same  accuracy  bounds  at  sufficiently  high  signal-to-noise  ratios.  The  best  array  is  then  the  one 
with  the  lowest  threshold  signal-to-noise  ratio. 

The  arrays  to  be  compared  all  have  seven  elements.  The  first  is  a  circular  array  with  elements 
evenly  spaced  in  angle.  For  a  CR  factor  of  1  wavelength,  this  array  has  a  radius  of  0.535  wavelength. 

The  second  array  is  a  spiral  array.  Its  elements  are  on  the  same  radials  as  the  circular  array, 
but  their  distance  from  the  center  grows  geometrically  by  the  factor  b.  This  parameter  will  be 
called  the  growth  factor  of  the  array.  In  this  example,  the  growth  factor  is  1.1.  For  a  CR  factor  of 
1  wavelength,  the  first  element  is  at  a  distance  of  0.391  wavelength  from  the  center. 

The  third  array  has  the  shape  of  the  letter  L.  There  is  an  element  at  the  intersection  and  three 
along  each  arm.  Along  each  of  the  two  arms,  the  element  positions  are  geometric,  i.e.,  r,br,b^r. 
For  radial  arrays  such  as  this,  the  parameter  b  will  be  called  the  growth  factor.  In  this  example, 
b  is  1.7.  For  a  CR  factor  of  1  wavelength,  the  distance  r  is  0.840  wavelength.  Although  several 
growth  factors  were  tried  for  the  spiral  and  L  arrays,  no  further  attempt  was  made  to  optimize 
these  parameters. 

The  three  array  geometries  are  shown  in  Figure  3. 

The  CR  factor  of  an  array  is  increased  by  expanding  the  array  or,  equivalently,  increasing 
the  frequency  of  operation.  Bounds  on  DF  error  for  the  circular,  spiral,  and  L-shaped  arrays, 
respectively,  are  shown  as  a  function  of  element  signal-to-noise  ratio  in  Figures  4  to  6  for  various 
CR  factors.  The  threshold  ESNR®  for  each  array  is  plotted  as  a  function  of  CR  factor  in  Figure  7. 

The  circular  array  clearly  has  the  lowest  threshold  ESNR,  while  the  spiral  is  somewhat  better 
than  the  L  for  most  CR  factors.  This  Figure  also  illustrates  the  general  principle  that  the  threshold 
ESNR  is  a  nondecreasing  function  of  CR  factor.  The  reason  for  this  is  that  increasing  the  size  of 
an  array,  and  hence  its  CR  factor,  causes  addition^  sidelobes  to  be  drawn  into  visible  space. 

It  must  be  emphasized  that  the  purpose  of  this  section  is  not  to  advocate  the  circular,  or  any 
other,  array  geometry  as  being  generally  superior.  It  is  simply  to  compare  three  particular  arrays 
on  an  equal  basis  using  the  WW  bound  as  a  tool. 


^Threshold  ESNR  is  defined  computationally  as  that  value  of  ESNR  at  which  the  multitest-point 
WW  bound  has  diverged  from  the  two-test-point  version  by  0.25  dB.  It  is  computed  to  a  precision 
of  0.2  dB  in  ESNR. 
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Figure  3.  Circular,  spiral,  and  L-shaped  arrays  having  CR  factors  of  unity. 
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Figure  4-  Bounds  on  DF  error  for  circular  array  for  various  CR  factors. 
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Figure  5.  Bounds  on  DF  error  for  spiral  array  for  various  CR  factors. 


26 


•  •  •  • 


•  • 


Figure  6.  Bounds  on  OF  error  for  L-shaped  array  for  various  CR  factors. 
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Figure  7.  Threshold  ESNR  behavior  of  the  three  arrays  vs.  CR  factor. 


12.  EXTENSION  TO  SWITCHED  SUBARRAYS 


One  possible  solution  to  the  problem  of  doing  reliable  DF  with  a  small  number  of  receiver 
chamnels  is  to  switch  the  available  channels  among  a  larger  number  of  antenna  elements.  The 
database  from  which  the  DF  is  derived  then  consists  of  several  looks  taken  from  each  of  several 
subarrays.  The  various  subarray  patterns  have  high  sidelobes  in  different  directions,  and  these  tend 
to  “average  out,”  producing  a  pattern  not  unlike  that  of  the  full  array. 

Calculation  of  the  WW  bound  for  switched  subarrays  requires  only  minimal  changes  in  the 
previous  analysis.  In  Equation  (9),  the  sum  that  appears  in  the  exponential  is  replaced  by  a  sum 
of  sums,  one  for  each  subarray,  and  the  first  factor  becomes  a  product  of  different  factors.  The 
covariance  matrix  C  is  different  for  each  subarray. 

The  effect  of  this  change  on  the  previous  analysis  is  simply  to  replace  Equation  (18)  with  a 
product  of  terms,  each  having  the  same  form  as  Equation  (18). 

5 

/lA/(h.,hj)  =  n[(A/./’+  l)detR,]-^*, 

5=1 


where  is  the  number  of  elements,  R,  is  the  matrix  of  Equation  (15),  and  N,  is  the  number  of 
looks,  all  for  the  sth  subarray. 

» 

The  test  points  in  the  switched  subarray  case  are  chosen  to  be  the  highest  sidelobes  of  the 
switched  subarray  pattern,  which  is  a  weighted  sum  of  the  subariay  patterns.  The  mathematical 
definition  of  the  switched  subarray  pattern  and  the  rationale  behind  this  definition  are  given  in 
Appendix  E. 

Figure  8  presents  the  results  of  a  Monte  Carlo  simulation  that  validates  the  WW  bound  I 

calculation  for  switched  subarrays.  The  simulated  array  is  a  7-element  spiral  with  a  growth  factor 
of  1.1  and  a  CR  factor  of  4.  Three  5-element  subarrays  are  used,  each  providing  eight  looks.  The 
three  subarrays  are  comprised  of  elements  (1,2,4,6,7),  (1,3,4,5,7),  and  (2,3,4,5,6),  element  1  being 
innermost.  As  before,  each  simulation  point  represents  an  average  of  1000  Monte  Carlo  trials.  At 
each  trial,  the  true  direction  cosines  Ux,Uy  are  chosen  at  random  from  the  unit  disk.  The  calculated  I 

bound  again  is  seen  to  give  a  reasonably  accurate  prediction  of  the  threshold  ESNR. 

Figure  9  shows  the  predicted  DF  accuracy  of  the  switched  subarrays  for  various  CR  factors 
(the  CR  factor  is  that  of  the  full  array).  This  can  be  compared  directly  with  the  unswitched  array 
accuracy  shown  in  Figure  5.  In  all  cases,  there  is  a  loss  of  about  2  dB  in  accuracy  at  high  ESNR 
(because  only  16  looks  are  available  from  most  elements  instead  of  24),  and  the  threshold  ESNR  is  • 

increased  by  approximately  2  dB.  Similar  results  are  obtained  for  the  circular  and  L-shaped  arrays. 
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AREA  OF  ERROR  ELLIPSE  (dB) 


Figure  8. 
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Comparisan  of  simulation  results  for  switched  spiral  array  with  WW  bound. 
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AREA  OF  ERROR  ELLIPSE  (dB) 


ESNR  (dB) 

Figure  9.  DF  error  bounds  for  switched  spiral  array. 
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13.  SUMMARY 


A  method  has  been  developed  for  comparing  sparse  array  configurations  of  ideal  elements 
based  on  DF  performance.  Arrays  to  be  compared  are  sized  so  that  each  has  the  same  DF  error 
bound  at  high  signal- to-noise  ratio.  They  are  then  ranked  in  terms  of  their  threshold  signal-to-noise 
ratios,  the  lowest  being  deemed  the  best.  The  method  makes  use  of  the  WW  bound  as  it  applies 
to  the  array  DF  problem.  This  bound  accurately  predicts  the  threshold  signal-to-noise  ratio  of  an 
ideal  array,  without  the  need  for  computer  simulation. 

As  an  example  of  its  use,  three  array  configurations,  a  circle,  a  spiral,  and  an  L,  are  compared. 

The  WW  bound  can  also  be  computed  for  a  slightly  more  general  situation  in  which  data 
is  available  sequentially  in  time  from  several  subarrays.  In  this  way  it  is  possible  to  collect  data 
from  a  large  array  with  a  smaller  number  of  multiplexed  receiver  channels.  As  an  example,  it  is 
shown  that  the  spiral  array  already  examined  can  be  operated  in  this  mode  with  only  a  small  DF 
performance  degradation.  Similar  results  are  obtained  with  the  circular  and  L  arrays. 

It  is  almost  certainly  true  that  the  threshold  ESN  R  of  a  real  array  is  determined  by  calibration 
errors,  rather  than  by  the  noise  and  small  sample  effects  that  the  bound  takes  into  account.  One 
should  not  expect  a  real  array  to  achieve  the  degree  of  DF  accuracy  predicted  by  the  bound. 
However,  it  is  hoped  that  if  one  array  is  better  than  another  in  the  sense  of  this  report,  and  they 
are  implemented  with  equal  attention  to  calibration  issues,  the  ordering  will  be  preserved. 


APPENDIX  A 

PROOF  THAT  EACH  ADDITIONAL  PERTURBATION  VECTOR 

INCREASES  THE  BOUND 


The  following  easily  verified  factorization  of  a  partitioned  matrix  is  valid  whenever  A  is 
invertible. 


A  B 
C  D 


I  0  \  /  A 


I  A-'B 


CA-‘  I  /  \  0  D-CA'*B  /  \  0  I 


Assuming  that  D  —  CA  ^B  is  also  invertible,  the  inverse  of  this  factorization  is 


A  B 
C  D 


I  -A-‘B  \  /  A 


0  I 


0  (D-CA->B)-‘  /  \  -CA-'  I 


I  0 


When  a  perturbation  vector  is  aaded,  the  H  matrix  becomes  (H  h),  while  F  becomes 


r  g 

a 


Use  of  the  above  factorization  for  this  matrix  gives  for  the  augmented  bound  matrix 


n'  =  (H  h) 


T 

g  a 


(H  h)^ 


=  (H  v)  (  ^  °  _  j 

\  0  (a-g^r-'g)  '  / 


where 


=  HF-'H^  +  cw^ 


V  =  h  -  HF-'g, 


C  =  {a-  g^r  ‘g)  ’ 

is  a  nonnegative  number  because  the  expanded  F  matrix  must  be  nonnegative  definite.  It  is  now 
evident  that 


x^n'x  >  x^Ox 


for  any  vector  x,  and  so  fl'  >  fi. 


I 


» 


I 


> 
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APPENDIX  B 

ON  THE  CONVERGENCE  OF  THE  WW  BOUND  TO  THE  RANDOM  • 

PARAMETER  CR  BOUND 


As  mentioned  in  Section  8,  it  has  been  claimed  that  if  only  incremental  test  points  are  used, 
the  WW  bound  becomes  the  random  parameter  CR  bound  in  the  bmit  as  t  — >  0.  The  example 
given  here  shows  that  this  is  not  always  the  case. 

For  simplicity,  attention  will  be  restricted  here  to  the  one-parameter  case,  so  that  the  bound 
is  a  scalar  quantity.  In  this  case,  for  a  single  incremental  test  point, 


2  _ f,u)}p _ 

E{[L^I'^{z\u  -i-  c,  u)  -  u  -  f,  u)]^} 

[/  dz  J  du  -I- f)p'/^(z,u)]^ 

f  dz  j  du  [p‘/^(z,  u  -f-  € )  —  p'/^(z,  u  -  €)Y 


(B.l) 


It  must  be  remembered  that  L,  defined  by  Equation  (2),  is  defined  only  for  values  of  u  where 
p{z,u)  >  0  and  so  the  u  integration  is  limited  to  this  region. 

The  integral  in  the  numerator  approaches  unity  as  t  goes  to  zero  (see.  e.g.,  Rudin  [11,  Thm 
9.5]).  In  the  denominator,  if  the  interchange  of  '  he  u  integration  and  the  limit  operation  can  be 
justified,  then  clearly  the  bound  approaches 


A  J  dz  j  du  (^p'/^(z,n))^  =  J  dz  J  du  (— In  p(z,  u))^p(z,  u) 


(B.2) 


which  is  the  CR  bound  for  random  parameter  estimation  as  given  by  Van  Trees  [1], 

On  the  other  hand,  suppose  that  for  all  z,p(z,  w)  is  nonzero  only  on  a  finite  interval  [a  <  u  <  6] 
and  is  nonzero  and  continuous  from  the  left  or  right  as  appropriate  at  the  end  points.  Then  the 
interchange  referred  to  above  cannot  be  Justified,  as  the  following  argument  shows.  The  u-integrals 
are  then  over  the  finite  interval  [a, 6],  and,  for  c  >  0,  the  denominator  of  Equation  (B.l)  is 

^  j  dz^j  du  p{z,u  +  ()  +  dup(z,u-f) 

-I--4  /dz  /  dw  [p‘^^(z,w -h  f)  -  p'^^(z,u  -  c)]^ 

^  J  Ja-^e 

Because  p(z,u)  is  not  zero  at  its  end  points  and  is  continuous  from  the  right  or  left,  as  appropriate, 
the  expression  within  the  large  brackets  is  of  order  £,  and  so  the  limit  of  the  first  term  above  as  c  |  0 
is  -1-00.  The  remaining  term  is  positive,  so  the  denominator  is  +oo  and  the  bound  goes  to  zero.  On 


» 


liJ 


•  •  • 
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the  other  hand,  Equation  (B.2)  may  well  give  a  finite  result,  as  for  example  if  z  is  a  complex  scalar 
and  p(2,ti)  is  defined  on  C*  x  (0, 1]  by 


p{z,u)  =  p{z\u)p{u)  =  -exp{-(z  - 

fi  z 

This  problem  can  be  avoided  by  choosing  the  a  priori  densities  so  that  p{z,u)  and  its  u- 
derivative  are  zero  at  the  ends  of  the  interval;  however,  this  may  greatly  complicate  the  numerical 
calculation  of  the  bound. 


» 


» 
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APPENDIX  C 

EVALUATION  OF  A  DETERMINANT 


» 


* 


•J 

I 

♦ 

The  determinant  of 

Q  =  I-c(xx^  +  yy^),  I 

where  c  is  real,  can  be  calculated  by  doing  an  eigendecomposition.  Because  Q  is  Hermitian,  it  has 
exactly  M  eigenvectors.  Af  —  2  of  these  are  orthogonal  to  the  subspace  spanned  by  the  vectors 
x,y  and  have  eigenvalue  1.  The  remaining  two  must  be  linear  combinations  of  x,y.^  To  determine 
these,  solve  I 


(I  -  cxx^  -  cyy")(x  +  oy)  =  A(x  +  ay), 
which  reduces  to 

(1  -  A  -  c(x^x  +  ax"y))x+  ((1  -  A)a  -  c(y"x  +  ay^y)y)  =  0. 

The  coefficients  of  x  and  y  must  both  be  zero.  The  value  of  a  that  makes  the  coefficient  of  y  zero 
is 


_ cy^x 

^  1  -  A  -  cy^y ' 

Substitution  of  this  result  into  the  remaining  coefficient  yields  a  quadratic  equation  for  1  —  A, 
(1  -  A)^  -  c(x^x  +  y^y)(l  -  A)  +  c^(x^xy^y  -  |x'^y|)^)  =  0, 


with  the  solution 

A  =  1  -  ^  ^x^x  +  y^y  ±  yj{x^x  -  y^y)^  +  4|x^yP  ^  . 
The  value  of  the  determinant  of  Q  is  the  product  of  the  eigenvalues. 


» 


i 


^If  X  and  y  are  colinear,  the  problem  reduces  to  a  simpler  one. 
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In  the  problem  addressed  by  this  report,  the  vectors  x,y  have  equal  length  so  that  the 
previous  equation  reduces  to 

A  =  1  -  c(M  ±  |x"y|), 

and  the  determinant  is 


det  Q  =  (1  —  Mcf  —  (c|x"y|) 


APPENDIX  D 
CALCULATION  OF  fi(h.,hj) 


The  calculation  of  5(h,,hj)  involves  evaluating  the  integral  in  Equation  (19)  with  p(u)  given 
by  Equation  (6).  The  region  of  integration  f/  is  a  disk  of  unit  radius  centered  at  the  origin,  and  the 
factors  in  the  integrand  are  nonzero  only  on  disks  of  unit  radius  centered  at  -h,,  — hj,  respectively. 
The  fundamental  problem,  then,  is  to  find  the  area  of  intersection  of  three  unit  disks. 

If  two  disks  of  unit  radius  intersect,  the  area  of  intersection  is  the  union  of  two  circular 
segments  (see  Figure  D-2).  The  inclusion  of  a  third  disk  of  unit  radius  can  have  four  possible 
results: 

(a)  Same  area  of  intersection  as  above  [Figure  n-2(a)j. 

(b)  Reduced  area  of  intersection  still  composed  of  two  circular  segments  [Figure  D-2(b)]. 

(c)  Reduced  area  of  intersection  bounded  by  thrte  circular  arcs,  which  is  the  union  of 
three  circular  segments  and  a  triangle  [Figure  D-2(c)]. 

(d)  No  intersection. 

The  area  of  a  circular  segment  (see  Figure  D-3)  that  subtends  angle  6  is 


Aseg  - 


^  -  sin  ^ 


The  angle  ff  can  be  expressed  in  terms  of  the  coordinates  of  the  corners  of  the  circular  segment  as 


0  =  2arcsin(d/2),  d=  J(x2-  Xi)2  +  (yj  -  yi)^. 


The  area  of  a  triangle  with  vertices  at  (zi.yi), (X2,y2),(a:3,y3)  is 


Atrt  =  ridet 


X2  -  Xi  X3  -  X] 

1/2  - 1/1  ys  -  yi 


These  formulas  make  it  possible  to  evaluate  Equation  ( 19)  without  resorting  to  numerical  integra¬ 
tion. 
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APPENDIX  E 

MAXIMUM  LIKELIHOOD  ESTIMATOR  FOR  SINGLE-SIGNAL 
DIRECTION-FINDING  WITH  NORMAL  AND  SWITCHED  ARRAYS 

Let  Zs,i,Zj,2, •  •  •  be  independent  A/,-dimensional  vector  measurements  taken  from  the 

sth  subarray,  s  =  1, . . . ,  5.  It  is  assumed  that  these  measurements  are  of  the  form 


where 


the  same  for  all  subarrays. 


=  Im, 


vf  (u)v,(u)  =  M„ 

for  all  u,  where  M,  is  the  number  of  elements  in  the  sth  subarray.  Assuming  that  the  condi¬ 
tional  probability  density  p(z|u)  of  the  measurements  is  zero  mean  Gaussian,  the  logarithm  of  the 
likelihood  function  for  the  measurements  taken  from  the  sth  subarray  is  [see  Equations  (8)  and  (9)] 

lnp(z,,i,z,,2, . . .  ,z,,N.|u)  =  -JV,[M,ln IT -f  ln(l -t- M,P)] 


l-l-Af.P 


v,(u)v"(u))z,,„. 


The  maximum  likelihood  estimate  of  u  is  that  vector  u,  which  maximizes  the  above  expression.  The 
only  term  that  depends  on  u  is 


'•<">= ttm:? 


the  expected  value  of  which  is 


E{r,(u)} 


— ^(]V,P|v"(uK(uo)P  +  M,). 

1  +  M,P 


(E.2) 


The  maximum  value  of  Equation  (E.2)  occurs  at  u  =  Uo-  The  array  power  pattern  5,(u,  Uq)  of  the 
sth  subarray  is  defined  to  be  the  dot  product  in  Equation  (E.2)  normalized  so  that  its  maximum 
value  is  unity,  i.e., 


ff,(u,uo)  =  |^v^'(u)v,(uo)|^ 

The  mean  value  of  the  log  likelihood  function  for  the  entire  set  of  switched  measurements  is  just  the 
sum  of  those  for  the  individual  subarrays,  each  of  which  is  meiximum  when  u  =  uq.  The  quantity 
analogous  to  Equation  (E.l)  for  switched  subarrays  is  therefore 


s 

r(u)  =  51! 

S=l 


(E.3) 


with  expected  value 


«!/ 


♦ 


The  switched  subarray  power  pattern  is  therefore  defined  to  be  a  weighted  sum  of  individual  subarray 
patterns,  the  weights  depending  on  the  strength  of  the  signal  as  well  as  the  number  of  elements 
and  looks  contributed  by  each  subarray. 


where 


s 

flf(u,Uo)  = 

5=1 


A, 


N,Mj  ^  N,MI 
1  +  M,P'  ^^l  +  M,P' 


I 
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For  simplicity,  attention  will  be  confined  here  to  subarrays  of  equal  size  (M,  =  M),  so  that 

5 

A,  = 

S=1 

i.e.,  the  weight  of  each  subarray  is  just  the  fraction  of  the  total  number  of  looks  it  contributes.  The 
sidelobe  locations  of  the  switched  array  pattern  are  used  as  test  points  when  evaluating  the  WW 
bound  on  switched  subarray  DF  error. 

The  estimate  u(z)  that  maximizes  Equation  (E.3)  is  the  maximum  likelihood  estimate  of 
direction  for  a  switched  subarray  DF  system.  It  is  this  estimator  that  was  used  in  both  the 
unswitched  and  switched  simulations  described  in  this  report. 


APPENDIX  F 

DERIVATION  OF  THE  INFORMATION  MATRIX  FOR  DF  WITH 

PLANAR  ARRAYS 


This  appendix  derives  an  expression  for  the  Fisher  information  matrix  F  for  the  array  DF 
problem  in  the  general  case  where  each  source  is  characterized  by  D  real  directional  parameters. 
The  directional  parameters  are  treated  as  constant,  unknown  parameters.  Although  £)  <  2  for 
most  DF  problems,  it  is  no  harder  to  handle  the  general  case,  and  the  results  may  be  useful  in 
other  applications. 

Assume  there  are  S  signals  (indexed  by  k,l)  each  having  D  parameters  (indexed  by  m,n). 
The  parameter  vector  for  signal  k  is  =  (“U.uzfci •  •  •  The  vector  u  is  the  concatenation  of 

these  vectors.  The  complex  5-vector  of  signal  amplitudes  is  assumed  to  be  a  zero  mean  Gaussian 
random  vector  with  known  covariance  matrix  P. 

The  Fisher  information  matrix  consists  of  5^  blocks,  each  of  dimension  D  x  D.  The  block 
corresponding  to  u/t,u;  is  (see  Van  Trees  [l]) 


dlnp(z\u)^T,d\npix\u)^ 


For  brevity,  denote  the  mnth  element  of  by  /„„.  Then  (see  DeLong  [12,  Appendix  B]) 

fmn  =  ^  Tr  (C-»  ^C-'  1^),  (F.2) 

du„M  duni 

where 

C(u)  =  I-t- V(u)PV"(u) 

is  the  multiple  signal  extension  of  Equation  (7),  in  which 
V(u)  =  (v(Ui),v(U2),...,v(us)) 

is  a  matrix,  the  columns  of  which  are  the  array  response  vectors  of  the  signals,  and  P  is  the 
covariance  of  the  Gaussian  signal  amplitudes.  The  necessary  derivative  is 


dUmk  OUrnk  dUmk 


Only  the  /:th  column  of  V(u)  depends  on  Umti  so  that 


ly  -  M  ^  I  ■  11 1  — 1^— 

«• 


aV(u) 


^V(Ufc)  Q 

’  ’■ 


..,0)  =  Vte^e*, 


(F.4) 


where  ^ 

V  =  ^  ^v(ufc)  ^(ut). 

*  duk  duik  '  du2k  dupk  ' 

and  Cm  and  are  unit  vectors  of  dimension  D  and  5,  respectively,  with  unit  element  in  row  m 
and  k,  respectively. 

Evaluation  of  Equation  (F.2)  using  Equations  (F.3)  and  (F.4)  gives 

/„„  =  iV  Tr  {C-«(V*e,„eJPV«  +  VPete3;V«)C->(V,e„e/'PV"  +  VPe,e^V")} 

=  iV  {(e[PV"C-'V,e„)(erPV"C-*V*e„  )  +  complex  conjugate  + 

(elVj^C-^  V,e„  )(ef  P  V"  C" '  VPe*  )  +  complex  conjugate  } 

=  2^»{(PW,)*„(PW*),„  +  (W*,)„„(PWP),*,}  (F.5) 

where 


w 

A 

v«c~‘v. 

(F.6) 

• 

W* 

A 

v^c-Wk, 

(F.7) 

W*, 

A 

Vj^C-'V,. 

(F.8) 

The  dot  notation  for  the  W  is  for  mnemonic  purposes  only;  it  does  not  indicate  differentiation.  * 

Let  (Aki)mn  and  denote  the  first  and  second  terms  inside  the  brackets  in  Equation  (F.5). 

Define 


W  =  (Wi,W2,...,W5) 


and 
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» 


•  •  •  •  • 


» 


^  W„ 

W,2  . 

■  W,s  ^ 

wi 

W2, 

W22  . 

•  Wjs 

^  W5, 

W52  . 

•  W55  y 

The  second  term  (Bfc/)mn  is  recognized  as  the  product  of  the  mnth  element  of  W*/  and  the 
scalar  c*;  =  ((PWP)^)*;,  so  the  contribution  of  this  term  to  the  complete  Fisher  information 
matrix  F  is 


B  = 


/  WiiCji  W12C12  ...  W15C15  ^ 

W21C21  W22C22  ...  W2SC25 


\  W51C51  Ws2C52 


WS5CSS  / 


=  w 


/  cii1£)1d 
C2i1£)1d 


f22lDlD 


C2S1d1d 


\  Csilplo  ‘^S21£)1d  •••  csslplp  / 


=  W  □  ((PWP)^®  Ipl^). 


(F.9) 


The  symbol  □  denotes  the  Hadamard  product;*  the  symbol  ®  denotes  the  Kronecker  product. 
The  first  term  can  be  written  as 


9 


(A*,)„n  =  (PW*)L(PW()*n  = 


®(AQB)mn  =  Omn^mn  The  Hadamard  product  is  defined  only  when  A  and  B  have  the  same 
dimensions. 

/  ojiB  ai2B  . . .  \ 


®A®B  = 


^21 B  ^22^ 


.  The  Kronecker  product  is  defined  for  any  two  matrices  A  and  B. 


\  . 


I 


» 
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showing  that  the  fc/th  block  of  A  is 


Ai/  =  (PWi)^e,efPW;,  (F.IO) 

which  is  a  dyad  involving  two  D-dimensional  vectors.  Any  two  such  vectors  x  and  y  satisfy  the 
identity 

xy^  =  {xl^)a(lDy^). 

Thus  Equation  (F.IO)  becomes 

A*,  =  ((PWfc)^e,lJ)0(lDerPW,). 


A. 


* 


A 


Define 

Qki  =  (PWt)^e,lS 

and  observe  that 

A*/  =  Q*;0(Q/ifc)^, 

so  that  the  complete  A  matrix  is  of  the  form 
A  =  QOQ^. 


The  full  Q  matrix  has  the  form 


'  (PW,)^e,lE  (PWO^ejl^ 

^  (PW5)^e,lS  (PW5)^e3lS 


(PWO^'e^l^  > 
(PW2)^e^lJ 


=  (  (PW)^e,lS,  fPWj^e^lS,...,  (PWl^e^lS) 
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•  •  •  •  •  •  • 


#. 


•  • 


=  (PW)^(eil]^,  CjlJ, 


The  matrix  on  the  right  above  can  be  written  more  compactly  as  Is0  Ip.  The  complete  A  matrix 
can  then  be  written  as 

A  =  ((PW)^(l5  ®  lS))0((l5  0  1d)PW).  (F.ll) 

Finally,  Equations  (F.ll)  and  (F.9)  in  combination  give  the  complete  Fisher  information  matrix 
for  the  directional  parameters, 


(F)uur  -  2N  3?{WD((PWP)^®1d1S)+((PW)^(I50T^)))D((I5®1d)PW)}.(F.12) 


If  the  array  elements  are  omnidirectional  and  only  one  signal  is  present.  Equation  (F.12) 
becomes  quite  simple.  P  and 

W  =  v"(I - - - vv")v  =  — — — 

^  A/P+l  ’  MP+l 

become  scalars,  and  because 

where  D  is  given  by  Equation  (22),  Equations  (F.7)  and  (F.8)  become 


W  =  >tv"(I  -  -^j^(^vv'')|vlD  = 


MP+l 


W  =  *^D^[v]"(I  -  ^^^w»)[v]D  =  *^D^(I  - 


The  first  matrix  inside  the  brackets  in  Equation  (F.12)  becomes 


k^MP^  T  P  T 

— iii_D^(I - - - 11^  )D, 

MP+l  '  MP+l  '  ’ 


while  the  second  becomes 
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» 


(MP+l) 


The  result  is  then  just 


F  = 


2Nk^MP'^ 

MP+  1 


Note  that  this  result  is  not  a  function  of  u. 


» 


» 


I 


» 


» 


» 


» 


» 


52 


i 


•  •  • 


•  •  • 


REFERENCES 


1.  H.L.  Van  Trees,  Detection,  Estimation,  and  Modulation  Theory,  Part  I,  New  York,  NY:  John 
Wiley  and  Sons  (1968). 

2.  E.  Weinstein  and  A.J.  Weiss,  “A  General  Class  of  Lower  Bounds  in  Parameter  Estimation,” 
IEEE  Transactions  on  Information  Theory  34,  338-342  (1988). 

3.  Ibid,  “Lower  Bounds  on  the  Mean  Square  Estimation  Error,”  in  Proc.  of  the  IEEE  73,  1433- 
1434  (1985). 

4.  A.J.  Weiss  and  E.  Weinstein,  “A  Lower  Bound  on  the  Mean-Square  Error  in  Random  Param¬ 
eter  Estimation,”  IEEE  Transactions  on  Information  Theory  31,  680-682  (1985). 

5.  D.H.  Johnson,  “The  Application  of  Spectral  Estimation  Methods  to  Bearing  Estimation  Prob¬ 
lems,”  in  Proc.  of  the  IEEE  70,  1018-1028  (1982). 

6.  H.  Cramer,  Mathematical  Methods  of  Statistics,  Princeton,  NJ:  Princeton  University  Press 
(1951). 

7.  D.F.  DeLong,  “Multiple  Signal  Direction  Finding  With  Thinned  Linear  Arrays,”  MIT  Lincoln 
Laboratory,  Lexington,  Mass,  Project  Rep.  TST-68  (13  April  1983). 

8.  F.M.  White,  “Performance  of  Bayes- Optimal  Angle-of- Arrival  Estimators,”  MIT  Lincoln  Lab¬ 
oratory,  Lexington,  Mass,  Technical  Rep.  654  (13  August  1984). 

9.  J.E.  Evans,  J.R.  Johnson,  and  D.F.  Sun,  “Application  of  Advanced  Signal  Processing  Tech¬ 
niques  to  Angle  of  Arrival  Estimation  in  ATC  Navigation  and  Surveillance  Systems,”  MIT 
Lincoln  Laboratory,  Lexington,  Mass,  Technical  Rep.  582  (23  June  1982). 

10.  R.J.  McAulay  and  L.P.  Seidman,  “A  Useful  Form  of  the  Barankin  Lower  Bound  and  Its 
Application  to  PPM  Threshold  Analysis,”  IEEE  Transactions  on  Information  Theory  15, 
273-279  (1969). 

11.  W.  Rudin,  Real  and  Complex  Analysis,  New  York,  NY:  McGraw-Hill  (1974). 

12.  D.  F.  DeLong,  “A  New  Algorithm  for  Adaptive  Radar  Signal  Processing,”  MIT  Lincoln  Lab¬ 
oratory,  Lexington,  MA,  Technical  Note  1979-72  (31  December  1979).  AD-A082856 


4 


i 


I 


9 


i 


REPORT  DOCUMENTATION  PAGE 


Foim  Apfimvtd 
0MB  Mo.  07O«-0fM 


10704^IMI.  WMmokn  DC  20000 


n  DOM  >  Sum  laM.  kmpim.  WA  oaoOAOQO.  and  ID  *•  OMOD  0 


■  M  OuOiaL  OWDI 


1  AGENCY  USE  ONLY  (Leavo  Wan/i[) 


2  REPORT  DATE 
16  Aunuot  1W3 


3  REPORT  TYPE  AND  DATES  COVERED 
Technical  Report 


4  TITLE  AND  SUBTITLE 

Use  of  the  Weiss- Weinstein  Bound  to  Compare  the  Direction -Finding  Performance 
of  Sparse  Arrays 


6.  AUTHOR(S) 


Darrol  F.  DeLong 


5  FUNDING  NUMBERS 


C  —  Fl%28-90-C-<)002 
PR  —  589 


7  PERFORMING  ORGANIZATION  NAME{S)  AND  AODRESSfES) 

Lincoln  Laboratory,  MIT 
P.O.  Bos  73 

Uxington,  MA  02173-9108 


B.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

TR-982 


9  SPONSORING/MONfrORING  AGENCY  NAME(S)  AND  ADORESS(ES) 

Department  of  Defense 

Ft.  George  C.  Meade,  MD  20755 


10.  SPONSORINGAXONfTORING 
AGENCY  REPORT  NUMBER 

ESC-TR-93-244 


1 1 .  SUPPLEMENTARY  NOTES 


12a.  DISTRIBUTION/AVAILABILITY  STATEMENT 


12b  DISTRIBUTION  CODE 


Approved  for  public  release;  distribution  is  unlimited. 


13.  ABSTRACT  (Mawman  200  tvo/nls) 

Weiss  and  Weinstein  have  recently  derived  a  new  lower  bound  on  random  parameter  estimation  error.  This  report 
makes  use  of  this  hound  to  develop  a  criterion  for  comparing  two  sparse  array  configyirations.  The  arrays  are  first 
adjusted  in  size  so  that  their  direction-finding  (DF)  accuracies  at  high  element  signal-to-noise  ratio  (ESNR)  are  the 
same  when  averaged  over  all  possible  directions  of  arrival  in  the  hemisphere.  Then  their  threshold  ESNRs  are 
calculated  using  the  Weiss-Weinstein  hound.  The  array  having  the  lower  threshold  ESNR  is  judged  the  better  array. 
As  an  example  of  the  appUcation  of  the  bound,  three  T-element  arrays,  a  circle,  a  spiral,  and  an  L-shaped  array,  are 
compared. 

The  bound  is  also  extended  to  the  case  of  switched  subarrays,  where  data  are  collected  from  several  subarrays 
sequentially  in  time  and  processed  as  a  unit.  The  purpose  is  to  achieve  DF  accuracy  approaching  that  of  the  full  array 
while  reducing  the  number  of  signal  processing  channels.  An  example,  using  five  receivers  on  the  7-element  spiral 
array,  shows  that  this  is  indeed  possible. 


14.  SL©JECT  TERMS 

Cramer-Rao  factor  element  n^nl-to-noise  ratio 

Weia6*WeineleiD  bound  twitched  arrays 

direction  finding 

15  NUMBER  OF  PAGES 
66 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 
OF  REPORT 

Unclassified 

18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

Unclassified 

19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

Unclassified 

20.  LIMITATION  OF 
ABSTRACT 

Same  As  Report 

NSN  7540-01 -280-5500 


Standard  Form  298  (Rev.  2-69) 
Prescribed  by  AMSI  Std.  239-18 
298-102 


♦ 


•  • 


> 


I 


•  • 


•  • 


1 


